Large Scale Genome-Centric Metagenomic Data from the Gut Microbiome of Food-Producing Animals and Humans

The One Health concept is a global strategy to study the relationship between human and animal health and the transfer of pathogenic and non-pathogenic species between these systems. However, to the best of our knowledge, no data based on One Health genome-centric metagenomics are available in public repositories. Here, we present a dataset based on a pilot-study of 2,915 metagenome-assembled genomes (MAGs) of 107 samples from the human (N = 34), cattle (N = 28), swine (N = 15) and poultry (N = 30) gut microbiomes. Samples were collected from the five Brazilian geographical regions. Of the draft genomes, 1,273 were high-quality drafts (≥90% of completeness and ≤5% of contamination), and 1,642 were medium-quality drafts (≥50% of completeness and ≤10% of contamination). Taxonomic predictions were based on the alignment and concatenation of single-marker genes, and the most representative phyla were Bacteroidota, Firmicutes, and Proteobacteria. Many of these species represent potential pathogens that have already been described or potential new families, genera, and species with potential biotechnological applications. Analyses of this dataset will highlight discoveries about the ecology and functional role of pathogens and uncultivated Archaea and Bacteria from food-producing animals and humans. Furthermore, it also represents an opportunity to describe new species from underrepresented taxonomic groups.


Background & Summary
The use of metagenomic approaches has revolutionized clinical microbiology allowing simultaneous identification of all potential pathogens without the need for culture-based methods 1,2 . For example, real-time metagenomic outbreak surveillance has also been useful in the identification and tracking of unknown infections, as such Shiga-toxigenic Escherichia coli (STEC) O104:H4 in Germany 3 and SARS-CoV-2 coronavirus 4 . In the clinical context, metagenomics is also a powerful weapon in the fight against antibiotic resistance pathogens in humans and animals 5 . From the use of advanced methods based on de novo assembly of metagenomic sequences, several studies have reported the importance of the resistome (e.g., collection of antibiotic resistance genes) 6 . Improvement in the identification and quantification of antibiotic resistance genes from complete or near-complete genes makes the assembly approach useful for characterizing novel antibiotic resistance genes and/or comparing them with well-known genes 7 . On the other hand, it is also possible to establish the link between taxonomy and functional annotation using long-assembled sequences 8 , which can improve the characterization of antibiotic resistance genes and the identification of pathogens.
It is well known that environmental microbiomes are hotspots of antibiotic resistance genes and that these genes can be exchanged between environmental and host-associated microbiomes 6 or between host-and host-microbiomes 9 . The One Health concept is a global strategy to study the relationship between human, animal, and environmental health. The exchange of pathogenic and non-pathogenic microorganisms among these settings, associating the interconnection between humans, animals, and the environment, has been the main focus of one health study 10 . For example, Mosites and collaborators 11 reported that human and animal microbiomes share the same species of their gut microbiome in rural livestock-owning households in western Kenya. Another study, conducted by Sun et al. 11 , demonstrated that the three-month exposure of students to livestock farms resulted in high sharing of antibiotic resistance genes and the microbial community. However, to the best of our knowledge, no one health data based on large-scale sampling and high-throughput sequencing by focusing on microbial genome reconstruction from metagenome data has been available in public repositories.
Here, we present a large-scale genome-centric dataset based on a pilot-study of 2,915 metagenome-assembled genomes (MAGs) from 107 samples (Supplementary Table 4). Data can be reused to test new hypotheses about the potential exchange of microbes between food-producing animals and humans or explored in the biotechnology, evolutionary, functional, or ecological context.

Methods
Data generation. Data was generated from GUARANI (One Health Brazilian Group) network. Initially, the aims of the GUARANI network's project were to quantify the abundance and diversity of antibiotic resistance genes (e.g., resistome) of a large number of samples in Brazil (South American), distributed in the major five Brazilian geographical regions (North region -Castanhal, 1°17′46.3776″ S-47°55′8.6016″ W; South Region -Blumenau, 26°55′10″ S 49°3.967′ W; Southeast Region -Bragança, 22°57′9.7″ S-46°32.651′ W; Midwest Region -Dourados, 22°13′16″-S 54°48.334′ W; Northeast Region -Fortaleza, 3°43′2″ S-38°32.584′ W), and to investigate the relationship between human and food-producing animal microbiomes, and the potential exchange of pathogenic and non-pathogenic microbes between these systems (Fig. 1A) by metagenomic approaches. Supplementary Table 1 describes information about sex, species, age of animals, and demographic localization of the farms and cities where samples were collected. In general, the experimental design was based on general and descriptive traits. To cover a high number of samples from all geographical regions of Brazil and have great potential to perform large-scale genome-centric metagenomic data, we choose to use the same samples in Illumina high-throughput sequencing. For this, 107 samples [humans (N = 34), cattle (N = 28), swine (N = 15), and poultry (N = 30)] were collected in triplicate from farms located in the five Brazilian geographical regions (Fig. 1B). For each region, properties were selected based on the criterion of simultaneous swine, poultry, and cattle rearing. Human samples were collected from healthy individuals who lived in the closest urban areas to the rural properties. The World Health Organization defines health as "complete physical, mental, and social well-being", in this study, we followed this concept to define adults (>18y-o) without any physical disease or infirmity as healthy individuals. All human data was anonymized, and the authors affirm that human research participants provided informed consent for the publication of the microbiome data and all information was approved by the research ethics committees. Data collection was approved by the Research Ethics Committee (CEP), Committee on Ethics in the Use of Animals (CEUA) from Universidade Federal de São Paulo (UNIFESP) and National System of Genetic Resource Management and Associated Traditional Knowledge SISGEN (Process numbers: 3.116.383, 2607170119 and AA1668A, respectively). (CEP and SISGEN). All Cattle, swine, and poultry samples were collected only from adult animals. In the sample collection, a swab was introduced in the first 2 cm of the rectal region to collect faecal samples of animals. Invasive rectal swabs were used only to collect samples from animals (swine, cattle, and poultry). For humans, the subjects were instructed to collect stool samples using a sterile fecal collection container with no preservative. A sterile charcoal swab was introduced in the stool specimen, followed by the rapid removal of stool excess by pressuring the swab against the container wall. The samples were stored and shipped to a central lab for DNA extraction.
DNa extraction and sequencing. DNA extractions were carried out under sterile conditions in a microbiological vertical laminar airflow hood. We did not use negative control samples (e.g., "blank swab") because the reagent and laboratory contamination were most problematic in low microbial biomass microbiomes (e.g., placenta or lung human microbiome) compared that find in high microbial biomass microbiomes 12,13 , as such that found in the faecal samples used in this study. DNA was extracted directly from swabs using the ZymoBIOMICS (Zymo, USA) DNA Miniprep Kit. DNA integrity and quantification were performed using a Qubit ® 2.0 Fluorometer (Thermo Fisher Scientific, AU). All samples were quantified by Qubit and organized on the sequencing plates according to the DNA concentrations obtained (Supplementary Table 2). The samples that had the same range of amount (ng) of DNA were in the same plate, since the number of PCR cycles of amplification of the libraries depends on the amount of initial DNA, according to Illumina protocol. The samples from the different hosts were treated together with maximum attention to avoid cross contamination. In short, sequencing libraries were prepared with the Nextera DNA Flex Library Preparation Kit (Illumina, USA) according to the manufacturer's protocol. Sequencing was carried out in the NextSeq. 500 System (Illumina, USA) using NextSeq. 500/550 High Output Kit v2.5 (300 Cycles), generating 2 × 150 bp reads.
Metagenome assembly, binning, and genome quality control. To increase the throughput and maximize the number of MAGs in this dataset, we choose a strategy based on co-assembly. This strategy has been used in several studies, including in the reconstruction of genomes from poultry 15 , cattle 16 , and human 17 metagenomes. In this case, samples were merged using the combination of host and region samples (See Supplementary  Table 3 to check each Co-assembly dataset). Metagenomes were assembled using Megahit software 18 with the meta-large option (-min-count 2-k-list 27,37,47,57,67,77,87). A total of 4,861,910,960 high-quality reads were used to assemble 1,676,286 contigs greater than 2,500 bp ( Table 1). The binning approach was used to reconstruct genomes from metagenomes based on the compositional traits of individual contigs (e.g., tetra-nucleotide frequency and coverage) using Metabat2 with default parameters 19 . We considered only the genomes that passed rigorous quality control to remove spurious and contaminated genomes in the downstream analyses. Genomes with completeness ≥50.0 and contamination ≤10.0 were used in the downstream analyses, following the Minimum Information about a Metagenome-Assembled Genome (MIMAG) of bacteria and archaea standards 20 in CheckM software 21 with CheckM (lineage workflow). A total of 2,915 MAGs were reconstructed (Table 2 and  Supplementary Table 3). Of these MAGs, 1,273 are high-quality drafts (≥90% of completeness and ≤5% of contamination), and 1,642 are medium-quality drafts (≥50% of completeness and ≤10% of contamination) (Fig. 2). The mean and standard deviation of genome size were 3.1 ± 1.4 Mbp, while the number of contigs had a mean of 263 ± 263. In addition, the mean genome size is compatible with those described in human stool communities 22 . On the other hand, we assembled contigs greater than 2.02 Mbp in MAGs from poultry metagenomes, indicating the accuracy of the metagenome assembly. All MAGs were submitted under the NCBI database and post-processing through NCBI's Contamination Screen to remove adaptador and cross-species contamination. taxonomy prediction. We used standardized bacterial taxonomy based on genome phylogenomics proposed by Parks and collaborators 23 , using the GTDB-Tk v1.3.0 software 24 (classify_wf workflow) and the most recent version of the Genome Taxonomy Database (GTDB) Release 05-RS95 23 . This workflow has been used to infer the taxonomy of MAGs, once improved classification of new uncultivated lineages and standardized taxonomy ranks based on the phylogenetic information. The most representative phyla were Firmicutes, Bacteroidota, and Proteobacteria (Fig. 3A), which are extensively studied in host-associated microbiomes 25 . However, many of the MAGs described here are potential new genera or new families (Fig. 3B), highlighting new insights about the www.nature.com/scientificdata www.nature.com/scientificdata/ ecophysiology of these new taxonomic groups. Regarding shared species between the four microbial community hosts, 45 genera were shared among distinct hosts (Fig. 3C -Supplementary Table 5). This includes environmental species with ecological importance in the digestive microbiomes (e.g., Cellulomonas and Azospirillum). Furthermore, four shared genera were generically assigned as SZUA-444, SZUA-584, UBA1305, and UBA8346, demonstrating the importance of this dataset to explore new taxonomic groups.

technical Validation
Here, we reported 2,915 draft genomes assembled from host-associated metagenomes. Illumina metagenomic reads used to assemble MAGs went through multiple steps of rigorous quality control, which included removing low-quality reads and host-associated sequences. Only a small proportion of the reads (14.64 ± 11.19%) were removed during the quality control, which had 0.22 ± 2.12% of host-associated reads (Supplementary Table 2). In a total, 4,861,910,960 high-quality reads were used in the downstream analyses.
A total of 37,755,059 contigs were generated during the metagenome assembly steps, being 1,676,286 contigs greater than 2,500 bp were assembled (Table 1). Small contigs (≲2,500 bp) were discarded because they carried less compositional signatures (as such used in the binning step: tetranucleotide frequencies and coverage) and can bias the construction of clusters during the metagenome-assembled genomes reconstruction step 28 . The longest contigs showed a mean of 1,083,245 ± 295,772 bp (max: 2,020,273; min: 690,014), demonstrating the effectiveness of the high sequencing depths used here. These results are similar to those already described in other studies reconstructed contigs greater than 900,000 bp using host-associated microbiomes like rumen metagenomes 29 or caecum chicken microbiome 15 .
Each metagenome-assembled genome (MAG) was validated using the rigorous standards defined by the Minimum Information about a Metagenome-Assembled Genome (MIMAG) of bacteria and archaea consortium 20 , considering only medium and good quality genomes assigned by the number of single-copy genes within a phylogenetic lineage 21 . Furthermore, only 33 (1.13% of the total dataset) MAGs showed adaptor or cross-species contaminations during the NCBI's Contamination Screen, demonstrating the high quality of this  Table 1. Number of reads and metagenome assembly metrics of each individual data set. Table 2. Number and quality of metagenome-assembled genomes (MAGs) of each individual dataset. 1 Genomes with completeness => 50.00 and contamination = <10.00 2 Genomes with completeness = >50.00 and = <90.00 and contamination = <10.00; 3 Genomes with completeness = >90.00 and contamination = <5.00.